Supersolid of indirect excitons in electron-hole quantum Hall systems 
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We investigate the ground state of a balanced electron-hole system in the quantum Hall regime 
using mean- field theory and obtain a rich phase diagram as a function of interlayer distance d and the 
filling factor within a layer. We identify an excitonic condensate phase, an excitonic supersolid phase, 
as well as uncorrelated Wigner crystal states. We find that balanced electron-hole system exhibits 
a supersolid phase a wide range of filling factors, with different crystal structure ground states. We 
obtain the ground state stiffness in the excitonic phases and show that the phase transitions from a 
uniform condensate to a supersolid is accompanied by a marked change in the stiffness. Our results 
provide the first semi-quantitative determination excitonic supersolid phase diagram and properties. 

PACS numbers: 



Introduction: Over the past decade, Bose-Einstein con- 
densation of indirect excitons [l[ in optically pumped and 
doped electron-hole double quantum wells has been ex- 
tensively explored 0, H, 0, @, 0] The concurrent theo- 
retical analysis has largely focused on either a uniform 
condensate, where a macroscopic number of excitons oc- 
cupy the zero-momentum state, or the condensation in a 
trap where the translational symmetry is broken explic- 
itly 0,1, M,M EH. o n the other hand recent obser- 



vations of the supersolid phase in 4 He 12J, ll3| have re- 
vived the interest in and questions about excitonic Bose- 
Einstein condensates with spontaneously broken symme- 
tries and their properties [3, [lj| . ( "Supersolid" phase 
of cold atoms in optical lattices has been extensively dis- 
cussed, although, in that case, the translational symme- 
try is explicitly broken by the optical lattice.) 

Lozovik and Berman [l6[ first discussed the coherent 
charge-density-wave (CCDW) ground state of indirect 
excitons in electron-hole double quantum wells. Recently, 
based on general principles, it was shown that electron- 
hole systems must support a supersolid phase of exci- 
tons due to their dipolar repulsion 17J . Although a sim- 




ple qualitative analysis implies the existence of a super- 
solid phase in electron-hole bilayer systems, the quantita- 
tive determination of the phase-boundary and supersolid 
properties is made difficult by the dispersion of electron 
(hole) bands. 

In this paper, we study an electron-hole system in a 
strong magnetic field, where the kinetic energy of carri- 
ers is quenched [ijj]. We consider a system with equal 
electron and hole filling factors v e = Vh = vt/2 < 1/2. 
CCDW ground states of such a system have been inves- 




20 |- however, such states are destabilized by 
2lL I2H . We obtain the ground state phase 



tigated 
fluctuations 

diagram in the (oI,vt) plane (Fig. Q}. Particle-hole sym- 
metry in the lowest Landau level maps the electron-hole 
system at vt onto an electron-hole system at 2 — vt 
and implies that the phase diagram is symmetric around 
Vt = 1. We verify that this exact relation is satisfied by 
our mean-field results. The phase diagram of a closely re- 



FIG. f: (Color online) Ground state phase diagram of 
electron-hole quantum Hall system. For small d < d ci (re- 
gion I), the ground state is a uniform excitonic condensate; 
at large d > d C2 (region IV), the ground state is uncorre- 
lated triangular Wigner crystals. For d ci < d < d C2 the 
ground state is a supersolid of excitons, with triangular (re- 
gion II) and anisotropic (region III) lattice structures respec- 
tively. This state has spontaneous interlayer phase coherence 
as well as spontaneously broken translational symmetry. The 
inset shows that ground state lattice anisotropy j(vt) changes 
discontinuously at the boundary between regions II and III. 



lated system - bilayer quantum Hall system where carri- 
ers in both layers have the same polarity, near total filling 
factor i*r = 1 - has been extensively explored [n| . These 
systems have provided signatures of excitonic condensa- 
tion in interlayer tunneling and counterflow experiments 
at small d [IH and in the presence of a bias voltage [23[ . 

For an electron-hole quantum Hall system, we find that 
(Fig. [J): i) for small d < d Cl (region I), the ground state 
is a uniform excitonic condensate irrespective of the to- 
tal filling factor ii) for large d > d C2 (region IV), the 
ground state consists of uncorrelated triangular Wigner 
crystals in the electron layer and the hole layer iii) for 
d Cl < d < d C2 , the ground state is a supersolid of exci- 
tons, i.e. a state with spontaneous interlayer phase coher- 
ence and spontaneously broken translational symmetry. 
This supersolid has either a triangular lattice structure 
(region II) or an anisotropic lattice structure (region III) . 
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To study the robustness of phase coherence, we calculate 
the ground state stiffness p s {d,VT) by considering the 
mean-field energy dependence on an in-plane magnetic 
field. We find that the phase transitions with increasing 
d - from a uniform excitonic condensate to a supersolid 
to Wigner crystals - are accompanied by marked changes 
in the stiffness. 

In the following section, we briefly sketch the details of 
mean- field calculations. Then we discuss the phase dia- 
gram (Fig. |lj) focusing on the excitonic supersolid phases, 
and the dependence of stiffness p s (d) on the interlayer 
distance. We conclude the section with a comment on 
electron-hole quantum Hall systems at total filling fac- 
tors vt and 2 — vt. In the last section, we summarize 
our results. 

Microscopic Model: Let us consider a bilayer system with 
electrons as carriers in the top layer and holes in the bot- 
tom layer, in a uniform magnetic field B = B±z + B\\x. 
The magnetic field normal to the layers B± quenches 
the kinetic energy of the carriers, whereas the in-plane 
magnetic field Bu allows us to study mean- field states 
with a winding interlayer phase. We assume a vanish- 
ing interlayer tunneling amplitude, keeping in mind that 
the in-plane magnetic field cannot be gauged away for 
any nonzero interlayer tunneling. The Hamiltonian for 
the system in the lowest Landau level approximation is 
given by [l^l 



H 



n £ 



! (q)/ 5 ai CT i(q)/ 5 o 



(i) 



where a — e(h) denotes electron (hole), A is the area 
of the sample, T4e(q) = Vhh(<l) = 2ne 2 /(eq) is the re- 
pulsive intralayer Coulomb interaction (e ~ 10 is the 
dielectric constant of the semiconductor) and \4/i(q) = 
— V ee (q) exp(— qd) is the interlayer attractive Coulomb 
interaction. The momentum-space density operator in 
second quantization is 



1 



k%k 2 



(fci|e- iq - r |fc 2 ) 



(2) 



where c CT /. (c k ) is the annihilation (creation) operator 
for a particle in the lowest Landau level state \k) = \n = 
0, k), Nff, — A/(2nl%) is the degeneracy of a single Lan- 
dau level, and Is = \J hc/eB± is the magnetic length. 
We define the excitonic condensate operator as 



Peh(q) = jt- ^2(h\e 



- tq ' r |fc 2 )cc 



t 

h, — k 2 



(3) 



and recall that single-particle wave- functions for holes are 
obtained by complex-conjugation from the single-particle 
wavefunctions for electrons, (r|n, k)h = ( r \n,k)l- Fol- 
lowing standard procedure, we obtain the Hartree-Fock 



Hamiltonian 
Hhf 



d 



B 



crq 



Uaa(q)paa(q) + £/ CT <r(q)A™(q)] ■ (4) 



The first term, U aa (q) = [K(q) - V&(q)] p ro (-q) - 
V c {c\)pas (— q) contains the intralayer Hartree and ex- 
change, and interlayer Hartree contributions respectively 
(e = h), and second term C/ CTff (q) = — Vd(q)pcrs (— q) de- 
notes the excitonic condensate contribution. These di- 
mensionless contributions are given by V a (q) = l/qls, 
V b (q) = y/^/2I (q 2 l 2 B /4), VM) = -e-^Vaiq), and 



V d (q) 



d 2 P ;2 
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(5) 



To obtain the self-consistent density matrices, we in- 
troduce a two-component operator a\ = \c\ ■ , Ch.-k\ and 
define the 2x2 matrix Green's function l24| 



G(Q;r) 
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(6) 



The electron and hole density matrices (p aa {Q)) as well 
as the (complex) excitonic order parameter (p e h{Q,)} are 
then determined from the equal-time limit (r — * CP) of 
this Green's function matrix. In the Hartree-Fock ap- 
proximation, the equation of motion for G matrix in the 
frequency space is given by lil. 25. 2r3| 
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S e /i(Qi, Qj ) 
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where I is the identity matrix, Q^ = Q ± k^, keZ^ = 
ydB\\/B± represents the relative displacement of single- 
particle wavefunctions in the two layers due to the Lorenz 
drift caused by the in-plane field [26[ , and we have defined 
a modified Green's function 



Q{Q) 



G ee (Q) G e/l (Q+)e^ x «- £ / 2 



(8) 



The self-energies in Eq. ([7]) are given by S CTlC r 2 (Qj, Qj) = 
C/<T 1( r 2 (Qy)exp(iQ i xQjz;|/2) with Qij = Q^Qj. The 
modified Green's function can be expressed as 



G(Q;iuj n ) = ^ 



A fc (Q)Aj(0) ^ 

iu) n - Lu k 



(9) 



where X\(Q) = [V£(Q), U£(Q )] are the eigenvectors of 
the self-energy matrix with eigenvalue u>k, 

S 



3 



^ee(QijQj) P$ij 

S/, e (Q7",Qj) pSi 



eh(Qi, Qj ) 

Sfcfc(Qr,Q7) 



: WfeAfe(Qi) 



Afc(Qi) 
(10) 
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We notice that, in contrast to_a similar analysis for bi- 
layer quantum Hall systems 



251 . 1261 ] - the chemical po- r— . 



tcntial appears with opposite signs in the electron and 
hole self-energies. The Hartree-Fock ground state energy 
Ehf^b) is calculated using the self-consistent density 
matrix using Eq.Q. 

Equations I®-® provide the set of equations that are 
iteratively solved. We check that the resulting density 
matrix satisfies the sum-rule [251 
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(Q)| 2 + \Peh{Q)\ 2 } = Pee(0) = V T /2. (11) 



Results: In order to determine the phase diagram in the 
d — vt plane, we compare the energies of uniform exci- 
tonic condensate (p e /i 7^ 0, p oc <5 q ,o), Wigner crystal in 
each layer (p eh = 0,p ee = p hh oc 5^q), and supersolid 
(Peh 7^ 0, Pee = Phh oc Sq.o) in the absence of an in-plane 
magnetic field B\\ =0. We use simplified anisotropic lat- 
tice with two primitive lattice vectors 
citebrey2000 ai = (a, 6/2) and a 2 = (0,6), and define 
the lattice anisotropy 7 = b/a. The lattice constant a is 
determined by the constraint that the unit cell contains 
one electron-hole pair, and the optimal value of 7 is ob- 
tained by minimizing the mean-field energy. Note that 
the triangular lattice (7 = 2/y3 « 1.155) and a stripe 
lattice (7 — ► 0) are its special cases. (We have found that 
square lattices have higher energy in the supersolid and 
uncorrelated Wigner crystal phases.) 

Figure [1] shows the ground state phase diagram. At 
small values of d < d Cl , we find a uniform excitonic con- 
densate as the ground state over the entire range of total 
filling factor (region I). For such a state, the mean-field 
equations can be analytically solved and we obtain the 
layer densities p ee (0) = phh{0) = and the excitonic 
condensate order parameter p e h(0) — \/vt(2 — vr)/2. 
These analytical results satisfy Eq. (fTT|) and imply that 
the uniform phase coherence order parameter is the same 
for systems with filling factors vt and 2 — vt- The mono- 
tonic increase in d Cl with decreasing vt [HI is consis- 
tent with observed strengthening of excitonic condensate 
phase in bilayer quantum Hall systems with layer imbal- 
ance [23j. At large values of d > d C2 (which is the same 
as d ci for v T < 1/3), the ground state is uncorrelated 
triangular Wigner crystals in the electron layer and the 
hole layer (region IV). 

For vt > 1/3, we find that the ground state is a super- 
solid at intermediate values of distance d Cl < d < d C2 . In 
region II, the supersolid has a triangular lattice struc- 
ture (7 = 2/V3) whereas in region III, the lattice is 
anisotropic. (The anisotropic supersolid with small 7 is 
similar to the CCDW solutions discussed in Ref.[19j]). 
The boundary between two regions corresponds to a line 
of first-order phase transitions from a triangular lattice 
to an anisotropic lattice. The inset in Fig. Q] shows the 
ground state lattice anisotropy j(vt) along the line d Cl . 
For vt < 0.6, 7 = 2/y3 is constant, as expected for 
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FIG. 2: (Color online) Mean-field energy as a function of 
lattice anisotropy 7 for different values of cL/Ib at total filling 
factor vt = 0.62. For d/ls = 1.26, the optimal value of 
7 ~ 0.89 corresponds to an anisotropic supersolid, whereas 
increasing d/ls shifts the optimal value to 7 = 2/v / 3 = 1.155 
corresponding to a triangular supersolid. 



a triangular lattice; near vt ~ 0.6, 7 drops discontin- 
uously and then reduces monotonically with vt- This 
discontinuity in 7 indicates a first order transition from 
a triangular to an anisotropic supersolid. The reverse 
transition, from region III to region II, can be induced 
by increasing the interlayer distance d at a given value 
of vt ■ Figure [2] shows the mean- field energy as func- 
tion of lattice anisotropy 7 for increasing d. We see 
that the optimal value of 7 jumps discontinuously from 
7 w 0.9 (anisotropic lattice) to 7 = 2/V3 (triangular lat- 
tice). The triangular lattice structure in regions II and 
IV suggests that the phase transition between them is a 
continuous phase transition. On the other hand, differ- 
ent lattice structures in region III and IV imply that the 
transition between them will be a first order transition. 

An important property of a Bose-Einstein condensate, 
uniform or otherwise, is the energy cost associated with 
spatial variations of its phase. A uniform in-plane mag- 
netic field induces such a variation in a uniform excitonic 



condensate |27j. We calculate the stiffness, defined as the 
variation of ground state energy with an in-plane mag- 
netic field [28| 



As = t 



1 d 2 E HF {k B ) 
A 



dk% 



(12) 



from the ground state energy. Figure[3]shows dependence 
of the stiffness on interlayer distance p s (d) at vt = 0.4 
We see a change in the stiffness slope at d Cl , where the 
system changes from a uniform excitonic condensate to a 
triangular supersolid and d C2 , where it changes from the 
supersolid to uncorrelated Wigner crystals, along with 
the expected monotonic decrease with d. Since the trans- 
port properties (in superfluid or supersolid phase) depend 
on the stiffness, generically, we expect that the phase 
transitions will be visible in transport measurements. 

An electron-hole system at vt is mapped on to a sys- 
tem at 2 — vt after a particle-hole transformation on both 
layers. Figure 0] shows that the mean-field density ma- 
trices, obtained separately for vt = 0.5 and vt = 1.5, 
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FIG. 3: (Color online) Dependence of ground-state stiffness 
on interlayer distance pa{d) for total filling factors vt = 0.4. 
The slope discontinuities in the stiffness occur at the critical 
layer separations d ci and d C2 respectively. 



indeed satisfy /J ee (q, 2 - v T ) = <5 q ,o - Pee(q, v T ) and 
p e /j(q, 2 — vt) — Peh(<l, v t)- Therefore, it is enough to 
restrict ourselves to Vt < 1. 



existence of the supersolid phase in the intermediate dis- 
tance regime and the accompanying change in the phase 
stiffness - are generically valid. 

Our results present an example of a supersolid whose 
properties can be investigated starting from a micro- 
scopic Hamiltonian; in particular, the study of low- 
energy excitations and transport properties in the super- 
solid phase is possible. We point out that this excitonic 
system, although analytically tractable, is significantly 
different from the zero-field system in which, at low ex- 
citon densities, the excitons behave as non-interacting 
bosons. Nonetheless, an experimental verification (or fal- 
sification) of our predictions will provide a better under- 
standing of the supersolid phase and the properties of 
excitons in the quantum Hall regime. 
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FIG. 4: (Color online) Ground state electron density p ee (r) 
(top) and excitonic order parameter |p e h(r)| (bottom), mea- 
sured in units of l/(2irl%) at total filling factors vt = 0.5 
(left) and vt = 1.5 (right). The mean-field results satisfy 
Pee{r, v T ) = l-/9ee(r,2 — vt) and p eh (r, v T ) = p e h(r, 2- vt)- 



Discussion: The supersolid phase, with two sponta- 
neously broken continuous symmetries, has been a source 
of extensive investigations in 4 He 12J, [l3|, [lj| , but has not 
been experimentally realized in any other system. Our 
results predict that electron-hole quantum Hall systems 
exhibit a rich phase diagram including a supersolid phase 
robust over filling factors vt > 1/3. We show that the su- 
persolid phase exhibits a nonzero stiffness and obtain the 
dependence of stiffness p s (d) on interlayer distance. It is 
well known that mean-field approach overestimates the 
stability ordered states. However, it is not clear whether 
it overestimates the stability of one type of order (ex- 
citonic condensation) over another (broken translational 
symmetry); therefore, we believe that our conclusions - 
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